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This study examines the formation of a front over a 
mountain due to a pre-existing deformation field. The 
dependent variables are assumed to be independent of x. The 
hydrostatic Boussinesq equations are used with no diffusion 
of heat or momentum. A modified coordinate system similar 
to Phillips' sigma system is used. The model is bounded at 
the top by a rigid plane and periodic boundary conditions are 
used in the horizontal. The solutions show the formation of 
a front within a finite period of time that tilts toward the 
cold air. The frontogenesis experiments with a mountain of 
wavelength 3600 km showed negligible topographic effects. 
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INTRODUCTION 



I . 

Williams and Plotkin (1968) have shown that atmospheric 
frontogenesis occurs for large time limits using the quasi- 
geostrophic equation. Their results are unrealistic in the 
respect that the frontal zone does not tilt with height. A 
comparison of the linear and nonlinear hydrostatic primitive 
equations by Williams (1972) shows that the linear solutions 
are close to the solutions of the quasi-geostrophic equations. 
The nonlinear solutions are more realistic and show the forma- 
tion of frontal discontinuities in a finite period of time. 
Hoskins and Bretherton (1972) showed that dynamic processes 
in addition to horizontal deformation can produce a front. 
Bannon (1983) and (1984) investigated quasi-geostrophic and 
semi-geostrophic frontogenesis over topography. Monotonically 
sloping topography beneath a deformation leads to the forma- 
tion of a surface discontinuity in the absence of an initial 
horizontal thermal gradient. A front associated with a 
translating deformation field over an isolated orographic 
feature is affected by the topographically forced divergence 
field as the front passes over the moutnain. Convergence 
will increase frontogenesis and divergence will decrease it. 

The purpose of this study is to investigate the effect 
of topography on a front which arises from deformation. 
Periodic boundary conditions are introduced because they 
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simplify the mountain flow simulation. The deformation is 
independent of height. The Boussinesq equations are used 
in which the dependent variables are independent of x. A 
coordinate system similar to that of Phillips (1957) is used. 
The equations are solved numerically and output is shown in the 
form of vertical cross-sections in the y-£ plane. 
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II. BASIC EQUATIONS 



The basic equations used in this study are the Boussinesq 
equations which neglect the compressibility of the atmosphere. 
The hydrostatic approximation is made and heating and friction 
are neglected. These equations will be applied to a model 
in which the flow is bounded above by a rigid horizontal 
plane . 

The Boussinesq equations [Ogura and Phillips, 1962] may 
be expressed in the z-coordinate system as 




-V$ - f k xv 



( 1 ) 




0 



( 2 ) 




(3) 




0 



(4) 



where 



K 



6 




(5) 




( 6 ) 



and 



14 



(7) 



7T 



0 




p 



In this model, the Coriolis parameter is constant, and the 
only quantity which varies, with longitude is $ . 

To simplify the finite differencing of the equations, a 
modified coordinate system is introduced in which the lower 
boundary surface becomes a coordinate surface. This system 
permits the application of the lower boundary condition at 
the surface of the mountain. The coordinate system is similar 
to the sigma system introduced by Phillips (1957). 

Following Walton (1968) and DeBoer (1970), the vertical 
coordinate, zeta , is defined as 



where D and z are defined in Fig. 1 and Z = D -z . With 

s r s 

these definitions, £ = 0 at the lower boundary and C = 1 at 
the upper boundary. 

Transformation of equations (1) through (4) from the 
z-coordinate system into the ^-coordinate system yields 



z - z 



s 



( 8 ) 



Z 





(9) 




0 



( 10 ) 



9 $ 

3C 




(ID 
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Fig. 1. Schematic diagram of the zeta-coordinate system 



( 12 ) 



V- (Zv) + Z H- = 0 . 



In equations (9) through (12) , V is the horizontal gradient 
operator on zeta surfaces. 

It is desirable to put the equations in flux form in order 
to satisfy certain computational properties. Equations (9) 
through (12) become 





(13) 





(14) 




0 



(15) 




(16) 




(17) 



where M 



uZ , M 



vZ and H = 6Z. 



x 



y 



The boundary conditions are 



0 



at 



C = 0 



(18) 



and 



C 



0 at 



C = 1 • 



(19) 
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For this study, the flow will be independent of x and periodic 
in y. $ is the only variable allowed to vary in the x- 

g 

direction, so - — = 0 for all other variables. 

o X 

In order to obtain an expression for $, equation (16) 
may be integrated with respect to £ giving 



$ = / 0Z dc + C . ' (20) 

0 e o 

The constant of integration, C, can be eliminated by taking 
the vertical average of equation (20) and subtracting from 
equation (20) which gives the following 



$ = -2_ J 0 zdc - / 6 ZdC) + V . 

U 0 0 0 



( 21 ) 



In the above equation,' 'represents the vertical average of 
a given variable. 

A value for $ must be found before a solution for $ can 
be determined. First, take the vertical average of (17) and 
apply the boundary conditions (18) and (19) which gives 



3 HT 5 

7T— M 
3y y 



= - Z c 



= 0 



( 22 ) 



0 



This equation shows that M must be a constant in order to 
satisfy the boundary conditions. Next, take the vertical 
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average of equation (14) and differentiate with respect to 
y which yields 



atay ( V 



ay 



2^ 



= I— (- + TrS-CZe) (C-l)^' - 



3 Z 



ay 



ay 



o 



ay 



’fM’ ) 
x 



( 23 ) 



The first term in equation (23) is identically zero from (22) 
and the resulting equation becomes 



l? (z ff> ■ fy ( " l? (M y vi + ^ (ze,( ^ 1, l - f z "> ’ (24) 

This equation can be solved diagnostically for This equa- 

tion must be applied in order to satisfy the boundary condi- 
tions (18) and (19). 

9 

To obtain an expression for - — , which is a constant, 

9 x 

equation (21) is differentiated with respect to x which 
yields 



9$ 

3x 



_a_ 

e o 



z 



(f ^ 
J o ax 






r li 
' J ax 



d?) + 



3J> 

3x 



(25) 



Since 



39 

ax 



= 0 in this study. 



3j> _ _3j>' 

3x 3x 



( 26 ) 



and this quantity is assumed to be independent of y. 
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Consider the domain averages of equations (13) and (14) 



which are 



d_ 

dt 




_ 7 3^> 

Z 3x 



~ + f M 



( 27 ) 



it ( V 





(Z0) U-l) 



Fz 1 

3y 




(28) 



. . 3 0 

where represents the y average. Initially, is computed 

d X 

geostrophically as 



S'?' f rp 

3x ^ M y t=0 



(29) 



If this relationship is used for all time and if the mass 
flux in the x-direction is zero initially, then will be 
zero for all time. The first two terms on the right side 
of equation (28) represent the mountain drag. The sum of 
these two terms for mid-latitude flow has been shown to be 
negative in the mean. This will lead to a continual loss of 
y-momentum. This kind of loss would make a quasi-steady final 
state impossible. An artificial source of momentum is required 

to maintain the mean meridional flow. This artificial source 

n 

9 

of momentum is provided by holding constant as obtained 
in equation (29). In this case, if is decreased below 
the initial geostrophic value, M will increase and restore 
the momentum loss. Thus, in the long term mean, the mountain 
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terms are balanced by the mean Coriolis force. However, for 
the horizontal scales to be considered in this study, this 
momentum source may not be needed since synoptic scale moun- 
tains cause very little wave drag. 

The dependent variables are to be divided as follows: 



u (x, y , C , t) 


✓ — X 
>1 

>< 

D 

II 


+ u(y,c,t) , 


v(x,y,C,t) 


= V d (x,y) 


+ v' (y , C,t) , 


• 

C(YrCft) = 


c ' (y,c,t) 


/ 


e(y,c,t) = 


0' (y,c,t) 


/ 


<f>(x,y,c,t) 


= $(x,y) 


+ <t>' (y,c,t) . 



The basic wind deformation field is chosen to be periodic in 
y and to have zero vorticitv and divergence. Thus, the 
deformation field is given by 

U d (x,y) = -D d Lsinh(^) sin (^) , (31) 

V d (x,y) = -D d L cosh ( jr) cos (^) , (32) 

where D- and L are constants. This deformation field was 
d 

used by Stone (1966) in his study of quasi-geostrophic 
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f rontogenesis . Since the flow is assumed to be independent 
of x in this study, equations (31) and (32) are evaluated at 
x = 0 resulting in 



U d (0,y) = 0 , (33) 

V d (0,y) = -D d Lcos(^) . (34) 



Substituting the expressions in equations (30) through (34) 
into equations (13) through (15) and dropping the prime 
notation yields 



„ ~ 3(V,Z) . ~ 

z + v d + u 9 y + a7 (vZu> + ac (5uZ) 



= - Z |^ + f vZ , 

3x 



(35) 



rs r\ 3 (V, Z ) «\ r. 

Z H + l? (vZV) + v ~W~ + v d z l? + 3? (CUZ> 



= - z 



^ + #- 6 Z(c-l) |^ - fuZ , 



3y 6 



0 



3y 



(36) 



z tl + t^ze) + f^z) 



= 0 , 



(37) 



where = -D^ Lcos (^) . This last set of equations is the 
set that needs to be solved for u, v and 0. In this study, 
equation (24) does not have to be solved if 5P is restored 
to the initial constant value at every timestep. 
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III. FINITE DIFFERENCE FORMULATIONS 



Staggered grids are used in both the horizontal and 

vertical. The vertical grid is divided into K-l levels of 

constant A£ (Fig. 2) . Quantities defined for the layers will 

be denoted by integer subscripts and the £ coordinate value 

will be denoted by £1^.. The quantities 0, u, v and $ are all 

defined on the integer levels. The vertical motion, £, on 

the other hand, is defined on the top, bottom and interface 

levels. These levels will be denoted with half-integer 

subscripts and the vertical coordinate will be identified by 

£k + iy 2 . T ^ e horizontal grid is divided into J-l sections of 

constant resolution (Fig. 3) . Momentum variables are defined 

on the integer subscripted points denoted by Y ^ . The momentum 

variables include u, v, M , M and V.. The mass variables, 

' x y d 

0, 4>, Z and £, are defined on the half-integer subscripted 
points. These points will be denoted as Va lues 

needed off their respective grid points will be identified 
by a caret and will be defined as 



S . 
3 




i-l )S i-l/2 



( 38 ) 



for quantities defined on the half-integer points and as 



<K1 +1/2 



Oki b u , v . * 

<K1 - V2 ° < 



<k + 1/2 D < 

<U ° u, v e, <P 

<k-i/2 G < 



<3/2 °< 

<>1 □ U, V, d, <p 

< 1/2 ° < 



Fig. 2. The vertical grid, showing the indexing 

convention and staggering of the variables 
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>* 

X 







+ 

*0- N ^ a >- 



<N l> ° > 




*0- N ^ a >- 
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Fig. 3. The horizontal grid, showing the indexing convention 
and staggering of the variables 



s 



j + 1/2 




for quantities at the integer points. Periodic boundary 
conditions are used in this study for all quantities. 

Equations (11), (13) through (15) and (17) are solved by 

using finite differences in y, £ and t. The finite differen- 
cing scheme is chosen to conserve the sum of kinetic energy 
and potential energy. 

The time differencing scheme uses a leapfrog scheme with 
a periodic use of the Euler backward scheme. The Euler 
backward scheme is used to eliminate the computational mode 
arising with the leapfrog scheme [Haltiner and Williams, 

1980) . 

Equations (11), (17) and (13) through (15) may be written 

in finite difference form as 



M M ' 




/s 3u 




YI. , /O - YI. 1 tn 

3+1/2 3-1/2 



(CZ) 



j ,k+l/2 




2 ' s j ,k-l/2 2 



+ 



C k+l/2 C k-l/2 




( 40 ) 
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M jyi „ M ,M __ J.- 

3 V . i ( y j+l,k y j/k } ( j+l,k + j,k } _ ( y j,k y j-l,k } ( j,k j-l,k } 



j 9t 



1A 



YI . , , /0 - YI . , /0 

3+1/2 3-I/2 



V i k+l +V i k " v -i k +v i k _i 

'__2 /k+1/2 _ 2 [ 21 ^ j/k-1/2^ 2 [ 

^k+1/2 " C k-l/2 



/ + f-Csi) (c-l) k V/2-V/2 

1 YI j+l/2 3-1/2 °0 3 ' K 3,K Yi j+l/2 YI j-l/2 



fU j,k Z j ,k ' 



(41) 



90 ^ 

J j+1/2 9 t 



M 



z , _ + 



, 6 j+3/2,k +e j+1/2, k . _M 9 j+l/2,k +9 j-l/2,k 

V j+l,k 2 [ y j,k^ 2 } _ 

Y . , , - Y. 

3+1 1 



• , 9 j+1/2, k+l +6 j+1/2, k ,,-. 9 j+1/2, k +6 j+1/2, k-1 

j+1/2 .k+1/2 1 2 J j+1/2 ^-1/2^ 2 / 

C k+l/2 " C k-l/2 



= 0 , 



(42) 



^j+l/2,k+l ^ j+1/2 , k 



C k+1 r ’k 



-2. Z 

0 Q j+1/2 



9 j+1/2 ,k+l +G j+1/2 ,k 
2 



(43) 



and 



My -M 

j+1 /k y j,k _ 

Y . , _ Y . 

3+1 3 



-Z ( C i + 1 /2,k + l/2-^ + l/2,k-l/2 K (44) 

3+1/2 C k+l/2 C k-l/2 
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IV. RESULTS 



All numerical results to be shown use the following 
numerical values for the constants : 



a 


= 


12 K 


B 


= 


375 m 


D 


= 


9 km 


D d 


= 


10 ^ sec ^ 


f 


= 


-4 -1 

10 sec 


g 


= 


9.81 m/sec 


L 


= 


3600 km 


O 

CD 


= 


300 K 


36 


s 

3z 


= 


4 K/km 



The following finite-difference increments are used. 

At = 180 sec 

Ay = 4 6 km (46) 

AC = 1/27 

A. MOUNTAIN WITH NO DEFORMATION FLOW 

The initial conditions are obtained from the quasi- 
geostrophic equations for steady state flow over topography 
[Walton , 1968; DeBoer, 1970] . The initial potential 
temperature field is given by 
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80 



e (y / z , 0 ) 



0 + B 

s 



s sinh[a(z-D)] ,y, 2a ,y> 

8T sinh " (aD) " " COS ( L } COS ( L } 



(47) 



80 



The quantity — is the constant initial lapse rate and 

o Z 



a = 



eJ-(# 



80 1/2 



f L 0 0 8 z 



■) 



(48) 



Substituting equation (47) into the thermal wind equation 
yields 



8u 
8 z 



i <f> -¥ sin <f> 1 

f 0 q L 3 z s inn ( ot D ) L tt L 



(49) 



This equation can be integrated with respect to z to give 



u(y,z,0) 



0 O 



80 1/2 
c 



8 z 



cosh [a (z-D) ] y 

sinh(aD) 1/ 



- sin (f) (Z 

TT L 




(50) 



The initial v field can be shown to be 



v ( y , z , 0 ) = V + aVB 



cosh [a ( z-D) ] 
sinh (aD) 



cos (^) 



(51) 



The initial w and <p are obtained from equations (4) and (3) . 

It is necessary to rewrite the initial conditions in the 
zeta-coordinate system. From equation (6), equations (47), 
(50), and (51) can be rewritten as 
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6 (y,c / 0) 



o.o 30 s sinh [aZ (c-1) ] . 2a . ,y x /r . ox 

9 s + B T? - slnh(aD) Sln <L> Sln L ' (52 > 



u(y,c,0) 



= B(-g. !>) 1/2 °°qh.taz^-i)] sin ( y 



0Q 3Z 



sinh (aD) 



— sin (£) (c -i) D , 

TT L Z 



(53) 



and 



v (y , £ , 0 ) 



v + “VB S2f|igfl£ill cos® 



(54) 



These initial conditions are shown in Figs. 4-7. z, and <p are 
given by equations (10) and (9). 

The value of 0 g was computed using the formula 



0 



s 



9e s r D 
dz 2 



+ Z(c-l)]. 



(55) 



The profile of the mountain used in this study is given by 



z = B cos (^) . (56) 

S Lt 

and is shown in Fig. 8. 

The initial conditions show anticyclonic vorticity over 
the ridge and cyclonic vorticity over the valley. There is 
divergence on the ups lope and convergence on the downs lope. 
These perturbations die off with height above the surface. 
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Fig. 4. 6 (K) . Initial. Mountain 
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Fig. 5. u (m/sec). Initial. Mountain only 
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Fig. 6. v (m/sec) . Initial. Mountain only 
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Fig. 8. Profile of the mountain with maximum height 375 



Figs. 9-12 show the results at t = 72 hours. Since the 
original fields were analytic, only a slight adjustment to 
the mountain can be seen in the v, and hence, .the z, fields. 

B. DEFORMATION OVER FLAT TOPOGRAPHY 

For this case, the initial potential temperature field 
contains a horizontal variation. The initial field is given 
by 



9 0 

0 (y ,z ,0) = -^-(z " §) “ — cos (£) . (57) 

The initial u field is again derived using the thermal wind 
relationship and can be shown to be 

u(y,z,0) = - sin (^) (z - j) . (58) 

The initial v field is identically zero and thus, w is forced 
to be zero by equation (3) . 

It is again necessary to write equations (48) and (49) 
in zeta-coordinates as 

9 0 

9<y,C,0) = -jfU-UD - y°°s ( L> ' (59) 



and 



u(y,c,0) 



2a . ,y, , 

— sin (f ) ( c 

7T L 




(60) 
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Fig. 9. 0(K). 72 hours. Mountain only 
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Fig. 10. u (m/sec) . 72 hours. Mountain only 
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Fig. 11. v (m/sec) . 72 hours. Mountain only 
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Fig. 12. £ (x 10 sec ) . 72 hours. Mountain only 



v and C are zero initially and <f> is given by equation (9) . 
Figures 13 and 14 show the initial conditions and Fig. 15 
shows the deformation flow given by equation (34). 

The front is developed by 24 hours as shown in Figs. 16-19. 
The frontal discontinuity forms along the point where v = 0 
and tilts toward the colder air. Cyclonic shear can be seen 
along the frontal boundary. The relative vorticity is less 
than zero on the north side of the front and is greater than 
zero to the south. 

By 48 hours. Figs. 20-23 show the strengthening of the 
frontal discontinuity. The tilt remains the same as at the 
24 hour point. On the other side of the domain, the tempera- 
ture gradient is weakening, which indicates frontolysis 
occurring in that region. 

C. MOUNTAIN WITH NON-MOVING DEFORMATION 

The previous two cases are combined in order to investi- 
gate the combined effects of a deformation flow and a 
mountain. The initial fields are given by 



9 (y,C,0) 



0 + B 



s 



s sinh [gZ ( £-1 ) ] 
8 z sinh (aD) 



sin (^) — ~ cos ( , (61) 



u(y,£,0) 




cosh [aZ ( £-1 ) ] . ( y. 

sinh (aD) S1IU I/ 




( 62 ) 
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Fig. 13. 0 (K) . Initial. Deformation only 
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Fig. 14. u (m/sec). Initial. Deformation only 
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Deformation flow 
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0 ( K ) . 24 hours. Deformation only 
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Fig. 17. u (m/sec). 24 hours. Deformation only 
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Fig. 18. v (m/sec) . 24 hours. Deformation only 
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Fig. 20. 6 (K) . 48 hours. Deformation only 
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Fig. 21. u (m/sec). 48 hours. Deformation only 
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Fig. 22. v (m/sec). 48 hours. Deformation only 
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Fig. 23. £ (x 10 sec ). 48 hours. Deformation only 



v (y , £ , 0) = V + aVB 



cosh [aZ ( £-1 ) ] 
sinh (aD) 



cos 






(63) 



and £ is given by equation (10) . These fields can be seen 
in Figs. 24, 25, 26 and 27. The deformation is the same as 
in the previous case. 

As in the previous case, a tightening of the temperature 
gradient has occurred by the 24 hour point, although not as 
strong (Figs. 26-29). The front has moved farther to the 
north and the tilt is smaller than in Case B. Throughout 
the next 24 hours, the front continues to move to the north 
and passes over the mountain ridge (Figs. 30-33). The 
gradient weakens as the front passes over the ridge and the 
tilt is negligible. In this experiment, the front weakens 
because it is advected away from the f rontogenetic part of 
the deformation field. 

D. MOUNTAIN WITH MOVING DEFORMATION 

The next four cases to be considered start with the same 
initial conditions as in Case C (Figs. 24, 25, 26 and 27). 
However, in these cases, the deformation moves to the north 
at a speed of 10 m/sec. Each case examines a different 
starting place for the maximum northward deformation flow. 
The cases have the maximum flow over the valley, on the 
upslope side of the mountain, directly over the ridge, and 
on the downslope side of the mountain. 
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Fig. 24. 0(K). Initial. Deformation and mountain 
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Fig. 25. u (m/sec) . Initial. Deformation and mountain 
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Fig. 26. 0 ( K) . 24 hours. Deformation and mountain 
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(m/sec). 24 hours. Deformation and mountain 
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Fig. 28. v (m/sec) . 24 hours. Deformation and mountain 
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Fig. 29. : ('< 10 sec ). 24 hours. Deformation and mountain 
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Fig. 30. 0 (K) . 48 hours. Deformation and mountain 
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Fig. 31. u (m/sec). 48 hours. Deformation and mountain 
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Fia. 32. v (m/sec) . 4o hours. reformation and mountain 
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Fig. 33. C (x 10 sec ). 48 hours. Deformation and mountain 



1 . Maximum Deformation over the Valley 



The initial deformation is shown in Fig. 15. As in 
the previous cases, the front has formed by 24 hours (Figs. 
34-37) . There is a small tilt toward the colder air in the 
north and the front has moved farther up the mountain than 
in Case C. The upper level front has passed over the moun- 
tain ridge. 

As the deformation flow continues to move northward, 
so does the front. During the next 24 hours, the speed of 
the front slows down, but the gradient in the temperature 
field tightens (Figs. 38-41) and forms a discontinuity. The 
cyclonic shear around the front has increased and there is 
strong vertical motion in that region. 

A careful comparison of the potential temperatures 
at t = 24 hours between Figs. 24 and- 34 shows very little 
change between the flat experiment and the mountain experi- 
ment, except for the translation in Fig. 34. The same con- 
clusion holds at 48 hours. Bannon (1983,1984) has shown that 
mountain effects on f rontogenesis will be important when the 
mountain forced convergence is of the same order as the 
frontal convergence. Comparison of v in Figs. 6 and 18 shows 
that the mountain forced convergence is only about 25 percent 
of the mountain forced convergence at 24 hours. This ratio 
is even smaller later in the frontal evolution. With the 
large scale mountain used in this and later experiments, the 
effects of the topography are unimportant. 
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Fig. 34. 6(K). 24 hours. Maximum deformation over valley 
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Fig. 35. u (m/sec). 24 hours. Maximum deformation over valley 
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Fig. 37. C ( x 10 sec ). 24 hours. Maximum deformation over valley 
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Fig. 38. 9(K). 48 hours. Maximum deformation over valley 
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Fig. 39. u (m/sec). 48 hours. Maximum deformation over valley 



o 




UI3Z 



71 



Fig. 40. v (m/sec). 48 hours. Maximum deformation over valley 
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2 . Deformation Shifted One-Quarter Wavelength 



In this experiment, the deformation field- is shifted 
one-quarter of a wavelength to the north with respect to 
the initial fields (Fig. 42) . The topography is also unim- 
portant in this experiment. At t = 24 hours, there is only 
a slight indication of frontal formation in the temperature 
field, but a much stronger suggestion of the front in the u 
field. The u field, however, shows little tilt (Figs. 43- 
46) . Williams (1968) has shown that this phasing of the 
deformation field is inefficient for f rontogenesis . 

3 . Reversed Deformation Field 

As is shown in Fig. 47, the initial deformation flow 
is in opposite phase with the original deformation field. 

The front can be seen at 24 hours in Figs. 48-51. As with 
previous cases, the front tilts toward the cold air, but in 
this case, the tilt is toward the south. This is a result 
of the reversed deformation field. The fields are mirror 
images of the fields in the basic experiment except for the 
shift . 

4 . Deformation Reverse of Case 2 

The last case to be considered is the one in which 
the maximum deformation flow is the reverse of Case 2 (Fig. 

52) . The first indication of frontal activity can be seen 
24 hours later (Figs. 53-56) especially in the u field. 

Strong cyclonic activity can be seen here with no evidence 
of a tilt. The temperature field shows only a small tightening 
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Deformation when maximum is on the upslope side of the mountain 
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24 hours. Deformation shifted 1/4 wavelength 
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v (m/sec). 24 hours. Deformation shifted 1/4 wavelength 
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sec ). 24 hours. Deformation shifted 1/4 wavelength 
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Deformation flow when maximum is over the ridge 
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0(K). 24 hours. Reversed deformation field 
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Fig. 49. u (m/sec). 24 hours. Reversed deformation field 
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Fig. 50. v (m/sec). 24 hours. Reversed deformation field 
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Fig. 52. Deformation when maximum is on the downslope side of the mountain 
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Fig. 53. 0(K). 24 hours. Deformation reverse of case 
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Fig. 54. u (m/sec) . 24 hours. Deformation reverse of case 
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Fig. 55. v (m/sec). 24 hours. Deformation reverse of case 
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Fig. 56. c (x 10 sec ). 24 hours. Deformation reverse of case 



of the gradient from the initial conditions. This is very 
similar to Case 2, but the maximum shear is at the lower 
boundary. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



In this study, a numerical model of frontogenesis was 
constructed which included surface topography. Periodic 
boundary conditions were included and a periodic deformation 
field was introduced. The model was successfully tested with 
frontogenesis with no topography and for stratified flow over 
a mountain with no deformation field. Combined experiments 
with topography and a moving deformation field were carried 
out. The experiments showed negligible effects on fronto- 
genesis by the topography. 

The next step in this research is to introduce smaller 
scale topography which should lead to a noticeable effect 
as suggested by the frontogenesis theories of Bannon (1983, 
1984). Other effects including vertical wind shear could 
be included later. 
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